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Abstract 

In a real system the heating is nonuniform and a second-order 
phase transition into a broken symmetry phase occurs by propagation 
of the temperature front. Two parameters, the cooling rate tq and 
the velocity vt of the transition front, determine the nucleation of 
topological defects. Depending on the relation of these parameters 
two regimes are found: in the regime of fast propagation defects are 
created according to the Zurek scenario for the homogeneous case, 
while in the slow propagation regime vortex formation is suppressed. 
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1. Introduction. 

A common issue of particular interest both in cosmology and in condensed 
matter physics is an estimation of the initial defect density produced during 
a phase transition into a broken symmetry state 0, 0- In cosmology this 
is an important issue so far as the later transitions are concerned — for 
example in extended models where there are stable defects produced at the 
electroweak transition. It is probably not of much practical importance for 
the GUT transition (if our ideas about scaling are correct) because in that 
case so much time has elapsed that no trace of the initial conditions remains 
[||, |]. The same thing is true no doubt for some condensed matter systems, 
but it would be good to identify some situations in which predictions about 
the initial density (not just the later scaling value) could be tested. We really 
need to find cases (i) in which the transition is second-order, (ii) in which it 
is of first order and goes by bubble nucleation, and (iii) in which it goes by 
the 'false vacuum' becoming totally unstable, which has been called spinodal 
decomposition. All these cases correspond to different scenarios of the phase 
ordering below the phase transition, and the creation and evolution of the 
defects prior to the final establishment of long-range order (phase coherence, 
or generally the coherence of the Goldstone variables). 

The great thing about superfluid helium-3 is that it allows such a wide 
range of possibilities. Most of these cases, or possibly all of them, may be 
represented in different regions of the parameter space: the normal 3 He to 
3 He-B or normal 3 He to 3 He-A transition is second order, while the A-to- 
B transition is first order and proceeds in different ways depending on the 
pressure (and magnetic field) and the extent of supercooling ||. Also of 
course, it adds another example to the list of possible systems, including 
nematics 0, helium-4 [[/J and high- and low-temperature superconductors, 
but an example which is in important ways more relevant since the fermion- 
boson interaction in 3 He shares many properties of the quantum field theory 
in high energy physics M . This is why it is so valuable as an analogue of the 
early universe. 

One of the most interesting experiments, which can shed light on the 
phase ordering is the 'mini big bang' produced by thermal neutrons || [K]]. 
There the exothermic nuclear reaction n+ 3 He — * p+fH + 0.76 MeV produces 
a local Big Bang in 3 He-B — a region of high temperature, T > T c , where the 
symmetry is restored. The subsequent cooling of this region back through 
the second-order superfluid transition results in the creation of a network of 
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vortex lines. These seeds of vortex lines are grown by the applied rotation 
and measured with the NMR technique @. 

Theoretically the formation of defects in a second-order transition into 
a broken symmetry state, discussed in 0,0], corresponds to a homogeneous 
phase transition. In real experiments the heating is never uniform. For 
example, in the neutron experiments the heating above T c by the neutron 
and the subsequent cooling occur locally in a region of about 10 /im size. This 
results in a temperature gradient and thus in a propagating boundary of the 
second-order transition. The temperature gradient is present also in other 
experiments where defect formation is influenced by the details of the cooling 
through T c ; see e.g. Jy] where the type of 3 He-A vortices created depends 
on the cooling rate. One might think that the phase of the order parameter 
would be determined by the ordered state behind the propagating boundary 
and that vortex formation would be suppressed in this geometry. On the 
other hand, in the limit of very rapid motion of the boundary one comes to 
the situation of an instantaneous transition and the Kibble-Zurek mechanism 
should be restored. Here we discuss the criterion on the propagation velocity 
of a second-order phase transition which separates the two regimes. We shall 
illustrate it with the example of the experimentally studied normal 3 He to 
3 He-A and normal 3 He to 3 He-B transitions. 



2. Phase ordering in a spatially homogeneous transition 

A nonequilibrium phase transition into the state in which the symmetry 
U(l) is broken leads to the formation of an infinite cluster of the topological 
defects - vortices or strings [|I|. As a result, when the temperature crosses T c 
two length scales appear. One of them is the conventional coherence length, 
which diverges at T c 

(For superfluid 3 He, where the Ginzburg-Landau theory is valid, one has 
v = 1/2). Another one is the mean distance £ between the vortices in 
the infinite cluster; it defines the scale within which the phase of the order 
parameter is correlated. This scale divereges with time when the vortex 
cluster decays. We are interested in the estimation of the initial density of 
the defects, initial, i-e., at the moment when these two scales become well 
defined. 
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According to Zurek |2[], initial is determined by the cooling rate tq, which 
characterizes the time dependence of the temperature in the vicinity of the 
phase transition: 

e(t) = l-^&«- ■ (2.2) 

J-c T Q 

Well defined vortices are formed at the time when the regions within the 
Ginzburg-Landau coherence length become causally connected. Causal con- 
nection is established by the propagation of the order parameter. In the case 
of superfluid 3 He the propagation velocity of the order parameter can be es- 
timated as the velocity of spin waves which just represent the propagating 
oscillations of some components of the order parameter. Thus one has that 
the corresponding velocity also depends on time and is given by 

c(f) « c^f) , (2.3) 

where Cq is of order of the Fermi velocity vp ~ 10 4 cm. 
At the moment t ~ t coh when 

c(t')dt' , (2.4) 

Jo 

the regions within the coherence length £,(t co h) are already connected, while 
outside they are causally disconnected. So the phase of the order parameter 
is well defined within the coherence length, but the phases in regions outside 
£ do not match each other. This corresponds to well defined vortices with 
the separation initial — £(tcoh)- The time t coh after the transition when this 
happens is 

tcoh « y/W^ , T « ^ « ~ 10~ 9 S , (2.5) 

and thus the initial separation between the vortices in the infinite cluster is 

initial = at C o h )^U-) 1 " ■ (2-6) 

TO 

For the more general case, when the phase transtion is not necessarily de- 
scribed by the Ginzburg-Landau theory, one has initial ~ £,o( T Q/ T o) a - 

The further development of the infinite cluster, which leads to its final 
elimination, has been the subject of the intensive investigations in phase 
ordering kinetics (see e.g. the review 
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In the neutron experiments in 3 He-B the estimation for tq ~ R^/D, where 
Rb is the size of the bubble and D is the diffusion constant near T c . This 
gives tq ~ 10 /is, and initial ~ 10~ 4 cm. In the typical A-phase experiments 
the cooling of the A-phase through the transition from the normal liquid N 
occurs during tq ~ 10 3 — 10 4 s [|TTJ, which gives initial ~ 10 Cm. 



3. The two regimes for a moving front 

In the presence of a temperature gradient two parameters characterize 
the noneqilibrium phase transition. In addition to tq, which determines the 
time scale of the temperature change, one has now the characteristic length 
scale of the temperature 

1 IVTI 



A T r 



(3.1) 



Combining with tq this gives the velocity of the propagating temperature 
front: 

u r « — , (3.2) 
T Q 

which is thus the velocity of the propagating second-order transition. 

The homogeneous result Eq. (2.6) is obtained when the velocity of the 
front is large compared to c(t co h), so the causality argument does work. This 
gives the estmation for the critical value of the velocity 

v Tc ^^^c (^ . (3.3) 

tcoh TQ 

If Vt < vtci the slowly moving front dictates the phase of the order parameter 
and formation of vortices should be significantly suppressed compared to 
the case of a rapidly moving front. In 3 He-B neutron experiments one has 
A ~ Ri, ~ 10 /jm (where Rb is the maximum radius of the bubble of normal 
fluid above T c ) and vt ~ 10 3 cm/s. This is comparable with cq(tq/tq) 1 I a ~ 
10 3 cm/s so vortex nucleation is not suppressed by the moving interface in 
this experiment. 

In the A-phase experiments the typical A ~ 10 2 cm and vt ~ Atq ~ 
10 -2 — 10 _1 cm/s, which is much less than c(t co h) ~ 10 cm/s. Thus at 
this low Vt Eq. (2.6) does not hold since the phase correlation across the 
front occurs faster than due to the temporal change of T. So the creation of 
vortices is markedly suppressed. However even in this regime formation of 
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defects has been experimentally observed, the formation of planar solitons in 
3 He- A @. 



4. Phase ordering behind a slowly moving front. 

One may expect that in the whole range of v t the initial density of vortices 
immediately behind the front is determined by the general scaling law: 

~^^t,t q ) =x a F{x , y) %= r_Q y = v_r 

where F(u) is some function, which has different asymptotes in the two 
regimes discussed above. From the previous sections it follows that for 3 He 
one has a = (3 = 1/4. The regime of fast propagation of the temperature 
front corresponds to the asymptote F(u) — > 1 when its argument u 3> 1- Let 
us find the asymptote in the regime of slow propagation, i.e. F(u <C 1). 

Since the propagation velocity of the order parameter c(t) slows down 
very near the transition, in some neighbourhood of the front, in a layer of 
thickness Az « t> T At, the situation becomes 'homogeneous'. Here At is 
obtained from 

r At (At) 3/2 
v T At^ I c{t)dt « c ^-^- , (4.2) 



which gives 



At m tq—j . (4.3) 

c o 

Outside this layer the condensate phase is already fixed due to the phase 
correlations with the low temperature regions, which are transferred by order 
parameter waves propagating along the vertical axis. So the only source of 
the mismatch of the phase originates within this thin layer, which means that 
the magnitude of the coherence length £(T) within the layer determines the 
initial distance between the vortices as a function of vt at slow transition: 

UW ~ at « At) « £o— , ^ < . (4.4) 

vt c t q 

Thus in the limit « C 1, we find F(u) ~ so the initial length scale is 
essentially larger than in the case of a rapidly propagating front, where: 

L^M « 6(^) 1/4 , ^>(^)^ 4 . (4.5) 

TO C Tq 
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Estimation of initial f° r a slowly propagating front in the A-phase gives 
6nitiai ~ 0.01 — 0.1 cm. This is significantly smaller than the size of the cell, 
~ 1 cm, and thus well explains the appearance of solitons during the cooling 
into the A-phase in [13| . 

On the other hand this is of order of the intervortex distance in the ro- 
tating cryostat and thus can influence the vortex texture which appears in 
the rotating cryostat when the superfluid transition occurs under rotation 



n]. In a field of 10 mT two types of vortices are competing in the equi- 
librium rotating state: singular one-quantum vortices with the core radius 
r CO re ~ £(T) and continuous two-quantum vortices (textures) with r core ~ £d, 
the dipole length £d ~ 10~ 3 cm. At low rotation velocities, Q < 2 rad/s (or 
in general at low vortex density), the array of singular one-quantum vor- 
tices has less energy. The experimental evidence is that well below 1 rad/s 
singular vortices are created after the superfluid transition under rotation, 
and no dependence on the cooling rate was observed. However, the exper- 
iment at higher velocity, Q ~ 1.4 rad/s, showed such dependence: at slow 
transition the equilibrium (singular) vortices dominate in the cell — their 
fraction is n s /(n s + 2n c ) ~ 0.8 — while at fast transition the fraction of 
continuous vortices 2n c / (n s + 2n c ) sharply increases from 0.2 up to 0.8. The 
change occurs in a jump-like manner at d{F ~ 6/iK/min which corresponds 
to tq ~ 3 • 10 4 s. One may expect that this is related to the phase ordering 
process, which leads to an initial vortex density larger than the equilibrium 
one. The further relaxation of the intial network towards equilibrium may 
discriminate between textures and singular vortices, which have significantly 
different scales. 



5. Conclusion. 

In a spatially inhomogeneous phase transition vortex formation depends 
on the velocity vt of the propagating front of the second-order phase transi- 
tion. There is a critical value of the front velocity, vt c ~ c o( t / t qY^ 4 : j which 
separates two regimes. For a rapidly propagating front, with vt > vt c , vor- 
tex formation is the same as in a spatially homogeneous phase transition. 
For a slowly propagating front, with vt < Vt C i vortex formation becomes less 
favorable whith decreasing vt- 
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